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A system of interacting bosons or fermions in a two-dimensional harmonic potential (or, equiva- 
lently, magnetic field) whose states are projected onto the lowest Landau level is considered. Generic 
expressions are derived for matrix elements of any interaction, in the basis of angular momentum 
eigenstates. For the fermion "ground state" (A'' = 1 Laughlin state), this makes it possible to exactly 
calculate its energy all the way up to the mesoscopic regime N ~ 1000. It is also shown that for 
TV = 3 and Coulomb interaction, several rational low-lying values of energy exist, for bosons and 
fermions alike. 



I. INTRODUCTION 

This paper is a sequel of Ref. where exact eigenstates were discussed for bosons with contact interaction in 
the lowest Landau level (LLL) of a strong magnetic field in two dimensions, as well as eigenenergies for fermions 
with Laplacian delta interaction, for which Laughlin wavefunctions are known to be exact eigenstates. In this paper, 
general expressions for matrix elements of an arbitrary central interaction — a sum of two body-interactions V{rij) 
whose Fourier transform admits a Laurent expansion v(k) = + ao + aik -\- . . . — projected onto the lowest Landau 
level are derived. These include 1/r" with any n, as well as contact (delta) and Laplacian delta interactions. An 
exact expression of the interaction energy E{N) for the A^-fermion "ground state" (the n = 1 Laughlin state, which 
is actually the ground state in the presence of a harmonic potential) is derived from which the large N asymptotic 
behavior can be obtained. For Coulomb interactions, the asymptotics is E{N) cx 7V^/^, which is confirmed by direct 
numerical calculation up to A^ = 1000. Also for Coulomb interactions, in the three-body problem, rational values of 
energy exist for low values of the total angular momentum, for bosons and fermions alike. 

Clearly, on the experimental side, we have in mind rotating Bose-Einstcin condensates [§| on the one hand, and 
strongly correlated Quantum Hall fermion droplets [3] on the other hand. In both cases a magnetic field is present, be 
it real in the quantum Hall case, or effective (due to the rotation of the condensate) in the BEC case. In the sequel, as 
a matter of simplification, we consider a harmonic trap one-body Hamiltonian, and the projection of the interaction 
is made on the one-body harmonic eigenstates 

(^,^10,0 =(^j ^'e-^— (1) 

Indeed, if a magnetic field were added to the harmonic trap, the one-body eigenstates corresponding to the LLL 
(Landau level number n = 0, angular momentum I > 0) would be the LLL-harmonic eigenstates basis (in complex 
coordinates) 

(z,z-|0,/)= ( z'e-^--" ^.,^0 (^^j ^'e-^— ^ (2) 
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where Wt = \/ '^'i + 1^^, i^c being half the cyclotron frequency. Note that since we diagonalize the system in a 
given angular momentum sector (the angular momentum operator commutes with the interaction Hamiltonian) , the 
magnetic field simply shifts the total energy by a constant term, which is therefore ignored here. 



II. MATRIX ELEMENTS 



The Hamiltonian for N interacting particles in a harmonic trap is 



i=l 



N 



N 



n = -2^didi + — ^ziZi+ vi\zi-zj\). 

i<j=l 



(3) 



i=l 



As long as the interaction potential V{\zi — Zj\) vanishes at infinity, so that the harmonic potential dominates, the 
asymptotics of the wave function can be detached as usually. 



tp{zi,zi,...,ZN,ZN) =exp ^-^ y^ZjZj^ 



x{zi,zi,...,zn,zn) ; 



then the Hamiltonian acting on x is 



where the free Hamiltonian Hq is 



N 



i<j=l 



N 



Hq = Yi-'^^i^i + + ^Zidi + OJZidi) . 



Prom now on one sets u) = 1. The LLL projector has the form 



Px 



^ n r 



X(-2l) ■2^1) •••) -ZjV) ^jv) • 



LLL functions are analytic 

for such a function, Px = X> and, since x does not depend on ^'s. 



N 



PHoX = uj{N + Y,^idi)X- 
Perform a Fourier transform of the interaction, 



and introduce the complex coordinates 



V{r) = / ^^^'Mk) ■ 



kr — I — , k — Kx -\- iKy 



One finally obtains the LLL-projected interaction 

N N 



i<j=l 



i<j=l 
N 



'k{z^-Zj) 



v{k) 



X(...,^;i + i-,...,2j -i-,...)e 2 dkdk 



(4) 

(5) 

(6) 

(7) 
(8) 
(9) 

(10) 
(11) 



EE fc^(-+-)+^.(fc)e-Vdfc ^ i ^ : „ d....... .l..n. x(^i,-,^^) 

i<j=l rai,n2=0 •^'^ 



'ni!n2!(ni +n2)!22("i+"2) dz'^'dz]-^' 



N oo 



(12) 
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where one has expanded every term containing the z's into Taylor series, substituted k = /cexp(i0), and integrated 
over (/). 

In the last expression of (|T2|) the expansion coefficient k'^"^-^v{k)e~^ dk is well defined as long as v{k) admits 
a Laurent expansion v{k) = X]m=-i ^mk™' such that 







fc2"+ifc™e--dfc = 2"+^r(n+l + -^) (13) 



An elementary LLL A'^-body wave function is 

N 



i=l 

for bosons (fermions), it has to be (anti)symmetrized. It is an eigenfunction of the total angular momentum, L = 
^ii ^"^^ thus an eigenfunction of iJo- It is obvious from (|12l) that ^'X]i<j=i ^(ki ~ -^^il) conserves the angular 
momentum (as any central interaction should do), therefore it is enough to diagonalize it in each sector of given L. 

The states with the lowest absolute value of angular momentum — for brevity, we will refer to them as "ground 
states" (which they are if there is a harmonic potential), both of bosons, 

X^{N)^l, L^O, (15) 

and of fermions, 

N 

X^oiN) = n (z. - Zj) , L = , (16) 

i<j=i 

are nondegenerate with respect to angular momentum, which implies that they are both eigenfunctions of 
-^^i<j=i ~ bosons, this can be seen directly by looking at Eq. where only the n = term 

survives so that 

P f2 ni^. - z,\) = I v{k)ke-'^dk (17) 

i<j=i 

and thus 



E{N) = ^(^^ j v{k)ke-^<lk 
(the diagonal matrix element of the interaction in the momentum representation.) 



III. FERMION GROUND STATE 

In the Fermi case it is also possible to solve the eigenvalue equation 

N 



P J2 Vi\^^ - ^MoiN) = EiN)x^,iN) , (18) 
i<j=i 



where PJ2i<j=i ^il^i - is given in ((T^ . 

The key observation is that, since said state Xoi^) is known to be an eigenstate of PJ2i<j=i ^il^i — it is not 
necessary to calculate the whole LHS of (IT51) . Being a Vandermonde determinant Xoi^) rewrites as 



N 



{N) = l[^l-' + -- - , (19) 



1=1 
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tN 



z\ ^ on 



where the omitted (A^! — 1) terms come from antisymmetrization. It follows that the coefficient in front of Hz^i -'i 
the LHS of Eq. ([T5| is necesseraly E{N). Return to the second line of Eq. (|12p and let xi^i, zn) be the monomial 



n 



N 



=1 ^; 



Then a term with given i and j in the sum on the RHS of that equation will be a sum of monomials in 



each of which only the powers of Zi and Zj are different from pi and pj, respectively (and the power of any zi with 
I is still pi). Moreover, the sum of all powers of z's, which is the total angular momentum, never changes. 

Hence, there are only two cases when -PX]i<j=i ^(l^i ~ ^j\)x can contain Hili ^i^^ o^'^ of its terms: (i) pi — l — l; 
(ii) pi = I — 1 {I hi)', Pi — j ~ ^] Pj — i ~ ^ (i-C-, Zi and Zj are interchanged; in the Vandermonde determinant, 
the corresponding monomial comes with a minus sign). Moreover, in each of these two cases, for given ni and n2, no 
more than a single term in the binomial expansion of {zi — Zj)"i+"^ will yield the desired contribution. In case (i), 



that term is z^ 
in case (ii), it is z, 

N 



(so that the powers of zi and Z2 stay unchanged after differentiation followed by multiplication); 

[so that the power of z^, which is j — 1 in x, becomes j — 1 — ni+ni+i — j = i — 1 
X] likewise for Zj]. The maximum possible values of ni and n2 are the powers of Zi and 



ni+i-j n2+j 
^3 



, respectively, in x- Taking this into account and gathering all the coefficients, we obtain 



E{N)= f{N,k)v{k)ke--dk 



(20) 



where 



N 



f{N,k)= 

i<]=l 



I 3 



3 i 



.ni— 712—0 



with 



(fc) 
(fc) 



(i - ni)n-,{j - n2)r 



ni— n2— 



2(ni+n2) 



(21) 



(ni!n2!)2 

^_2)"i+"2+i-« 
ni\n2\{ni + I - j)\(n2 + 3 - ^y■ V2 



2(ni+n2) 



(x)„ = ^^^"^ is the Pochhamnier symbol. 

The summation over ni and 112 can be performed explicitly, by noting that 

Cijnin2{k) = fiini{k)fjjn2{k) , 
dijnin2{k) ( l)"* /iirii (^)/iin2 (^) i 



where 



and that 



Hence, 



^-^"^^)^ n!(i+'^-j)! (^''"^"(^ 



^/yn(fc) — TT TTT-i' ^ 



n=0 



P 

(i- l)!"j-i I 4 



(22) 



(23) 



(24) 



(25) 



N 



f{N,k)^ ^ 
i<i=i 



I t) (t 



For example for the first values of N one has 

fc2 



/(2,fc) 



2 

fc6 



1 



/('I,'!) 



25A:« 
'6144 3072 



fc6 



23/s4 
lY 



^l-l 



(26) 



6fc2 + 6 



(27) 



5 



Substituting into Eq. (1201) gives, within the LLL-projection approximation, the energy of the A'^-fermion Vandermonde 
state Xo (-^) foi' ^^ly central pairwise interaction. 

Further simphfication is possible. A recurrency relation is 



fiN + 1, fc) - f{N, k) = L]v_i Lm (^^ + 



4 



(28) 



where in the first term, X)p=o ^p(^) ~ ^n^^i^) has been used. This is much more efficient than (j26p . as it requires to 
compute a single sum for each subsequent iV, instead of a double sum. 
One wants to find an expression for the chemical potential 



Using 



H{N) = E{N + 1)-E{N) 



1 r/" 

1,2 /.nJV + 1- 



(29) 



(30) 



one rewrites the second sum in Eq. as J2iLo ^''i\(N+i)\ — ^(~*' 2 — i + N, k'^/A)^, which, taking into account 
yields 



^l{N) = 2 / v{2yG) 



^ N+l — i 



Ax 



(31) 



where U{a,b,x) is the confluent hypergeometric function of the second kind. 
Note that if v{k) is a power, 

v{k) ^ck"" , 

the integration in Eq. ((20|) can be performed explicitly, using [3] 



(32) 



to! n! 



r(a) 3i^2(a, a — A, 1 + 7 + m; a — A — n, 1 + 7; —1) . (33) 



As a result. 



Miv) = c2'"+^r(- + i) 



(2)7v^i(- f )jv ^ fm 
{N-1) 



~-o)N „ /TO TO TO 

, 3^2 I — + 1, — + 1, iV + 1: — 
1)!7V! ■^^V2 '2 ' '2 



- 7V + 1, 2: -1 



(34) 



N 



AT (^■^ + * ~ A^ ^ l)iv ^ /"m m TO 
_ (-l)"-' (.^l),(^+lJ), 3^.(y + l.y+A--» + 2,iV + l;--. + 2,A.-. + 2;-l 



(when u(fc) is a Laurent series, /x(A'^) can be obtained as a corresponding sum over to). 

From ((3T|) one can, in the large A^ limit, obtain the asymptotics behavior of /i(A), at least in the case of Coulomb 
interaction to = —1. The first term ^^e~'^^ LN{x)L]^_-j^{x) can be simplified using the asymptotics 



(35) 



where Ja{x) is the Bessel function, and i/ = 4n + 2a + 2. The second term yields a sum of integrals which converges 
to a constant ~ 0.9. As a result. 



poo 

N e' 
Jq 



'Jo{Vi^N + 2)x)JiiVWx)— - 0.9 

x 

The energy is obtained by integrating the continuous version of (|29| . di?/dA^ — /i(A) 

E{N) ~ c (^-^N^^^ - 0.9N^ . 



c I -VAf-0.9 

TT 



(36) 



(37) 
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FIG. 1: The function E{N) (dots) versus Eq. ((57)) (continuous curve). 



The N^^^ scaling is easy to understand. The number of pairs grows as iV^, whereas the characteristic radius of the 
system in the ground state, which is the radius of the classical orbit with imax = iV — 1, grows as \/N — and the 
same should be true of the mean interparticle distance. 

As an illustration, we have obtained exact numerical results for E{N) with the Coulomb interaction. A convenient 
normalization is c = ■y/2/7r, which renders the results rational. Remarkably, exact results can be obtained for up to 
N ^ 1000, for which a "brute-force" calculation, involving A^! terms, would clearly be impossible^. Partial results 
are shown below: 



N 



E{N) 



1 

2 

87 

64 

5147 
2048 

514 095 



131 072 

1977 801361250 785 



140 737 488 355 328 

1 859 029 096 417 154 793 530 197 844 505 235 



40 564 819 207 303 340 847 894 502 572 032 



A plot of the (discretized) logarithmic derivative, 

log E{N) - log E{N - 1) 



D{N) = 



logAr-log(A- 1) 



(38) 



in Fig. 2 is clearly consistent with limjv-).oo D{N) = |, in accordance with Eq. ([57)) . 



IV. THE THREE-BODY COULOMB CASE 



Coming back to Eq. (|T2)) . for Coulomb interaction, one has 



N 



N 



P vn h (-l)"^r(ni + n2 + ^) {z, - z,r^+-- 9"^+"^ 



i<] = l 



i<j — l ni ,712—0 
N 



ni!n2!(ni +7^2)! 2"i+"^+5 dz'^'dz''^' 



= 1 n=0 ^ 



(39) 



This can be directly diagonalized in each sector with a given number of particles N and angular momentum L, with 
the basis formed by (anti)symmetrized functions of the form ([14p . The general structure of the spectrum is similar 
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FIG. 2: The dicretized logarithmic derivative D{N) as a function of (dots) versus dE{N) /dlog N , with E{N) given by 
Eq. (|37|) (continuous curve). 



to the delta interaction case: There are centcr-of-mass excitations, so that above each state with energy E there is 
a "tower" with energies E + n, n = 1,2, .. . . Only the "pure relative" eigenstates, devoid of these excitations, are 
of interest. We restrict ourselves to the 3-body problem. For the bosons with delta interaction, all the 3-body states 
turned out to have rational energies [1] (with a suitable choice of an overall factor c). With Coulomb interaction, 
though, irrational values start appearing rather low in the spectrum, for bosons and fermions like. Scaling away the 
overall irrationality, as before, by putting c = •\/2/7r, one finds all the eigenvalues of "pure relative" states, up to the 
appearance of irrationalities, are, for bosons: 



L 


E 





3 


2 


33 
16 


3 


51 
32 


4 


1755 


1024 


5 


3153 


2048 


6 


3(27749±7V766249) 


65536 



and for fermions: 



L 


E 


3 


87 

64 


5 


4881 
4096 


6 


4119 
4096 


7 

8 


140283 
131072 

63255 
65536 


9 


3(10025047±3Vl07141413705) 


33554432 



V. DISCUSSION 

The opportunity to calculate the energy of an eigenstate of N interacting two-dimensional bosons or fermions is 
certainly due to the fact that the LLL projection simplifies the situation. It reduces the dimension of the single- 
particle phase space [Bl and even more importantly, the "ground state" (Bose condensate for bosons, Laughlin state 
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for fermions) ends up being an eigenstate of the interacting Hamiltonian, which means that ah one has to compute 
for that state is a single matrix element. Nevertheless, even this simplified setup has a physical meaning, which makes 
our results applicable to real systems. 

The relevant case is when the LLL is separated by a gap from the rest of the spectrum. This happens when the 
whole system rotates with angular speed w, or if there is a strong magnetic field. But if the LLL is fiat (which happens 
if there is a magnetic field but no harmonic potential) , all the LLL TV-body states have the same degenerate energy. 
Our result for the ground" state, with the minimum L, is valid (it still does not mix with the other states), but not 
meaningful physically, as that state is not separated by an energy gap from states with higher values of L. This 
changes if a harmonic potential adds L(u)t — i^c) to the energy. One can then claim that if the interaction is weak 
enough compared to the gap, the exact energy of the A^-body ground state is known. 

One has to be careful, however, when taking the thermodynamic limit. As soon as N{LUt — Wc) becomes bigger than 
ujc, the energy of the lowest single-particle state in the first LL becomes smaller than that of the A^-th single-particle 
state in the LLL. Actually, the LLL projection approximation breaks as soon as N(ujt — (^c) ~ i^c- Therefore, for our 
result to be interesting, the oo limit has to be taken first, and then the thermodynamic limit N oo. 

Finally the same techniques could be applied to excited states with higher values of L. To do so one woud have to 
evaluate the matrix element {p'ip'2 ■ ■ -PNlPYliKj^i '^il^i - ^j\)\PiP2 ■ ■ -Pn), where \p1p2 ■ ■ -Pn) = n;Ii -^f' properly 
symmetrized or \p1p2 ■ ■ - Pn) — Y\ii:j^i{zi — zj) H/li properly antisymmetrized. 
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